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Abstract 

Recent developments in the modelling of the dark matter distribution in our Galaxy 
point out the necessity to consider some physical processes to satisfy observational data. 
In particular, models with adiabatic compression, which include the effect of the bary- 
onic gas in the halo, increase significantly the dark matter density in the central region 
of the Milky Way. On the other hand, the non-universality in scalar and gaugino sectors 
of supergravity models can also increase significantly the neutralino annihilation cross 
section. We show that the combination of both effects gives rise to a gamma-ray flux 
arising from the galactic center largely reachable by future experiments like GLAST. 
We also analyse in this framework the EGRET excess data above 1 GeV, as well as the 
recent data from CANGAROO and HESS. The analysis has been carried out imposing 
the most recent experimental constraints, such as the lower bound on the Higgs mass, 
the b — ► S7 branching ratio, and the muon g — 2. In addition, the recently improved 
upper bound on B(B S — > has also been taken into account. The astrophysical 

(WMAP) bounds on the dark matter density have also been imposed on the theoretical 
computation of the relic neutralino density through thermal production. 
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1 Introduction 



It is now well established that luminous matter makes up only a small fraction of the 
mass observed in Universe. A weakly interacting massive particle (WIMP) is one of 
the leading candidates for the "dark" component of the Universe. Indeed, WIMPs may 
be present in the right amount to explain the matter density observed in the analysis 
of galactic rotation curves [1], cluster of galaxies, and large scale flows [2], implying 
0.1 < tlDNih 2 < 0.3. Actually, the recent data obtained by the WMAP satellite [3] 
confirm that dark matter is present, and they lead to the value 0.094 < Q,DMh 2 < 0.129. 

One of the most promising methods for the indirect detection of WIMPs consists 
of detecting the gamma rays produced by their annihilations in the galactic halo [4]- 
[30]. The amount of gamma-ray fluxes observed will depend of course on the nature 
of WIMPs through their annihilation cross sections, but also on their density. This 
is the reason why the inner center of our galaxy (~ 100 pc from the center), where 
the dark matter density is large, is the main theater of dark matter searches through 
gamma-ray signatures. Atmospheric Cherenkov telescopes or space-based detectors 
are used for this search. For example, one of the space-based experiments, EGRET, 
detected a signal [31] that apparently is difficult to explain with the usual gamma-ray 
background. Also, the atmospheric telescope CANGAROO-II claims a significant de- 
tection of gamma rays from the galactic center region [32]. In addition, the atmospheric 
telescope HESS has also observed gamma-rays [33] in Sagittarius A* direction. Pro- 
jected experiments might clarify the situation. This is the case of the space telescope 
GLAST [34], which is scheduled for launch in 2006, and will be able to provide a larger 
sensitivity. 

Concerning the nature of WIMPs, the best motivated candidate is the lightest neu- 
tralino, a particle predicted by the supersymmetric (SUSY) extension of the standard 
model [35]. In particular, in most of the parameter space of the minimal supersymmet- 
ric standard model (MSSM) the lightest supersymmetric particle (LSP) is the lightest 
neutralino. Thus it is absolutely stable and therefore a candidate for dark matter. 
In addition, the neutralino is an electrically neutral particle, and this is welcome since 
otherwise it would bind to nuclei and would be excluded as a candidate for dark matter 
from unsuccessful searches for exotic heavy isotopes. 

Recently, SUSY dark matter candidates have been studied in the context of realistic 
halo models including baryonic matter [36]. Indeed, since the total mass of the inner 
galaxy is dominated by baryons, the dark matter distribution is likely to have been 
influenced by the baryonic potential. In particular, its density is increased, and as a 
consequence typical halo profiles such as Navarro, Frenk and White (NFW) [37] and 
Moore et al. [38] have a more singular behaviour near the galactic center. The conclu- 
sion of the work in Ref. [36] is that the gamma-ray flux produced by the annihilation of 
neutralinos in the galactic center is increased significantly, and is within the sensitivity 
of projected experiments, when density profiles with baryonic compression are taken 
into account. 

Nevertheless, this result has been obtained assuming particular values for the neu- 
tralino mass and annihilation cross sections <jj, and therefore a more detailed analysis 
computing explicitly <jj in the framework of SUSY must be carried out. It has also 
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been discussed recently that another way of increasing the gamma-ray flux produced by 
the annihilation of neutralinos is to consider the possibility of non-universality in the 
scalar and gaugino sectors of the MSSM [20, 21, 26, 28]. In particular, in the context of 
supergravity (SUGRA) the flux can be increased significantly. This is so when depar- 
tures from the minimal supergravity (mSUGRA) scenario, where the soft terms of the 
MSSM are assumed to be universal at the unification scale and radiative electroweak 
symmetry breaking is imposed, are taken into account. Using non-universal soft sca- 
lar and gaugino masses the annihilation cross section of neutralinos can be increased, 
producing a larger gamma-ray flux [26, 28]. 

The aim of this work is to discuss in detail both effects, halo models with baryonic 
compression and non-universal soft terms in SUGRA, on the gamma-ray flux. We will 
see that a significant enhancement is obtained, and regions of the SUGRA parameter 
space are compatible with the sensitivity of present and projected experiments. 

Let us remark that in addition to the astrophysical bounds on the dark matter 
density discussed above, the most recent experimental constraints are also considered 
in the analysis. In particular, we implement the lower bounds on the masses of SUSY 
particles and Higgs boson, as well as the experimental bounds on the branching ratio 
of the b — > «7 process and on a^ USY , for which the more stringent constraint from e + e~ 
disfavors important regions of the SUGRA parameter space (see e.g. Ref. [23, 26]). 
In addition, we have also taken into account the last data concerning the B s — > 
branching ratio. It is now known that the upper bound on this process constrains 
strongly the parameter space of non-universal SUGRA models, and has important 
implications for direct searches of dark matter [39, 40]. Thus an analysis in the case of 
indirect detection is necessary and has been carried out here. 

The paper is organized as follows. In Section 2 we will discuss in general the gamma- 
ray detection from dark matter annihilation, paying special attention to its halo model 
dependence. In particular, the modifications produced in the different profiles due to 
the effects of baryonic compression will be studied. The particle physics dependence 
will be reviewed in Section 3 for the case of neutralinos in SUGRA. We will study 
the general case when scalar and gaugino non-universalities are present. In addition, 
the most recent experimental and astrophysical constraints will be taken into account 
in the discussion. In Section 4 both effects, halo models with baryonic compression 
and non-universal soft terms, will be considered for the study of the gamma-ray flux. 
The theoretical predictions will be compared with the sensitivity of present and pro- 
jected experiments, such as EGRET, CANGAROO, HESS, and GLAST. Finally, the 
conclusions are left for Section 5. 

2 Astrophysical inputs and adiabatic compres- 
sion 

As discussed in the Introduction, we are interested in the annihilation of dark matter 
particles in the galactic center [5]. There are two possible types of gamma rays that 
can be produced by the annihilation. First, gamma-ray lines from processes XX ~^ 77 
[10] and XX ~^ 1% [H]- This signal would be very clear since the photons are basically 
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mono-energetic. Unfortunately, the neutralino does not couple directly to the photon, 
the Feynman diagrams are loop suppressed, and therefore the flux would be small [12]. 
For more recent analyses of this possibility sec Rcfs. [9, 17, 26]. 

On the other hand, continuum gamma rays produced by the decay of neutral pions 
generated in the cascading of annihilation products will give rise to larger fluxes [13]. 
We will concentrate on this possibility in the following. 

Let us then discuss the theoretical formulas for the computation of the continuum 
gamma rays, and in particular the astrophysical inputs one has to use. As we will 
see below, the gamma-ray flux can be significantly enhanced for specific dark matter 
density profiles. 



2.1 Gamma-ray flux 

For the continuum of gamma rays, the observed differential flux at the Earth coming 
from a direction forming an angle ip with respect to the galactic center is 



dN l 1 f 

0> 7 (^) = ^-^(a l v)-— I / p 2 dl , (2.1) 

dtiy 6irm x J line D f sight 



where the discrete sum is over all dark matter annihilation channels, dN^/dE^ is the 
differential gamma-ray yield, is the annihilation cross section averaged over its 

velocity distribution, m x is the mass of the dark matter particle, and p is the dark 
matter density. Assuming a spherical halo, one has p = p(r) with the galactocentric 
distance r 2 = I 2 + Rq — 21Rq cos tp, where Rq is the solar distance to the galactic center 
(~ 8 kpc). It is customary to separate in the above equation the particle physics part 
from the halo model dependence introducing the (dimensionless) quantity 

m = ( 7fTrv77 a V / P 2 {r(l^)) dl . (2.2) 

8.5 kpc V0.3 GeV/cmV June of sight 

Thus one can write 

$ 7 (£ 7 ,^) = 0.94 x 1(T 13 cm" 2 s" 1 GeV" 1 sr" 1 

" ?SGf^)(^)'™- w) 

Actually, when comparing to experimental data one must consider the integral of J(ip) 
over the spherical region of solid angle AO, given by the angular acceptance of the 
detector which is pointing towards the galactic center, i.e. the quantity J(ACl) with 



• ' I'^-hcldO. (2.4) 



must be used in the above equation. For example, for EGRET AO is about 10 3 sr 
whereas for GLAST, CANGAROO, and HESS it is 10~ 5 sr. 
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The gamma-ray flux can now be expressed as 

0.94 x KT 13 cm" 2 s" 1 

px dm ( ( ai v) \ A00GeV\ 2 =. . . ^ , n r . 

where Ethr is the lower threshold energy of the detector, Concerning the upper limit 
of the integral, note that neutralinos move at galactic velocity and therefore their 
annihilation occurs at rest. 



x 



2.2 Halo models 

A crucial ingredient for the calculation of J in (2.4), and therefore for the calculation of 
the flux of gamma rays, is the dark matter density profile of our galaxy. The different 
profiles that have been proposed in the literature can be parameterised as [41] 

where po is the local (solar neighborhood) halo density and a is a characteristic length. 
Although we will use po = 0.3 GeV/cm 3 throughout the paper, since this is just 
a scaling factor in the analysis, modifications to its value can be straightforwardly 
taken into account in the results. Highly cusped profiles are deduced from A-body 
simulations 1 . In particular, NFW [37] obtained a profile with a behaviour p(r) oc r~ l 
at small distances. A more singular behaviour, p(r) oc r~ L5 , was obtained by Moore 
et al. [38]. However, these predictions are valid only for halos without baryons. One 
can improve simulations in a more realistic way by taking into account the effect of 
the normal gas (baryons). This loses its energy through radiative processes falling to 
the central region of forming galaxy. As a consequence of this redistribution of mass, 
the resulting gravitational potential is deeper, and the dark matter must move closer 
to the center increasing its density. 

This increase in the dark matter density is often treated using adiabatic invariants. 
The present form of the adiabatic compression model was numerically and analytically 
studied by Blumental et al. [43]. This model assumes spherical symmetry, circular or- 
bit for the particles, and conservation of the angular momentum M{r)r = const., where 
M(r) is the total mass enclosed within radius r. The mass distributions in the initial 
and final configurations are therefore related by M,(rj)rj = [M^rj) + M£>M{ff)]ff, 
where Mj(r), Mt>(r) and Mdm(t) are the mass profile of the galactic halo before the 
cooling of the baryons (obtained through A-body simulations), the baryonic composi- 
tion of the Milky Way observed now, and the to be determined dark matter component 
of the halo today, respectively. This approximation was tested in numerical simula- 
tions [44, 45]. Nevertheless, a more precise approximation can be obtained including 
the possibility of elongated orbits [36]. In this case, the mass inside the orbit, M(r), 
is smaller than the real mass, the one the particle 'feels' during its revolution around 



1 For analytical derivations see e.g. the recent work [42], and references therein. 
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the galactic center. As a consequence, the modified compression model is based on the 
conservation of the product M(r)r, where r is the averaged radius of the orbit. The 
time average radii is given by x ~ 1.72x a82 /(l + 5x) ' 085 , with x = r/r s , and r s the 
characteristic radius of the assumed approximation. 

The models and constraints that we used in this work for the Milky Way can be 
found in Table I of Ref. [36]. We have fitted the resulting data of that work with 
the power-law parameterisation of eq. (2.6). The results are shown in Fig. 1, and 
summarized in Table 1. There we label the resulting NFW and Moore et al. profiles 
with adiabatic compression by NFW C and Moore c , respectively. As one can see, at small 
r the dark matter density profile following the adiabatic cooling of the baryonic fraction 
is a steep power law p oc r~ 7c with 7 C rs 1.45(1.65) for a NFW c (Moore c ) compressed 
model 2 . We observe for example that the effect of the adiabatic compression on a NFW 
profile is basically to transform it into a Moore et al. one (see Table 1). 

Let us remark that NFW and Moore et al. profiles can be considered as a lower and 
upper limit, respectively. For example, the one recently proposed by Diemand, Moore 
and Stadel in Ref. [46], where p{r) oc r~ L16 , is between both 'standard' profiles. It will 
be interesting for the analysis in Sect. 4 to use an average value of J for compressed 
scenarios. In particular, we will consider a value J' = 5 x Jnfw c > an d we will refer 
to this scenario as NFW^. for simplicity. It is worth noticing in this sense than some 
mechanisms can be advocated to reduce the effect of the compression like angular 
momentum transfer to dark matter or formation of the central black hole by spiraling 
and merging of two black holes. The study of such complex mechanisms is far beyond 
the scope of this work. 

Of course, these results have important implications for the computation of the 
gamma-ray fluxes from the galatic center. In particular, in Table 1 we see that for A$7 = 
10~ 3 sr one has </nfw c /^nfw — 145 and «/Moore c /^Moore et al. — 77. Thus the effect of 
the adiabatic compression is very strong, increasing the gamma-ray fluxes about two 
orders of magnitude. Similarly, for AQ = 10~ 5 sr one obtains </nfw c /^nfw — 953 and 

"^Moorec /^Moore et al. — 95. 

Let us finally remark that the J calculation has been regulated by assuming a 
constant density for r < 10~ 5 kpc. This procedure has consequences essentially for 
divergent J when 7 > 1.5. To check our results we have also calculated J in the fol- 
lowing way. We used a cut off in density, p(r < R s ) = 0, R s being the Schwarchild 
radius of a 3 x 1O 6 M black hole. In addition we also considered dark matter anni- 
hilation and estimate its effect by requiring an upper value for dark matter density, 
p(r < r max ) = p(r max ) = m x /((av).t BH ) ( t B H ~ 10 10 yr, being the age of the 
black hole). We also applied the same procedure with a more realistic cut-off value 
p(r < 10~ 6 pc) = suggested by dark matter particle scatterings on stars [47]. This 
kind of effects can be significant only for the Moore c halo model we consider in this pa- 
per and increase our J value by less than an order of magnitude for all the (av) values 
corresponding to our wide SUSY parameter range. In this sense our analysis is con- 
servative. More refine treatments of dark matter annihilation effects in the innermost 
region of the galaxy can be found in [48, 49]. 



2 It is worth noticing that we obtain the same order of magnitude as the authors of Ref. [48] . 
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Figure 1: Fit, interpolation and data for: a) 
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Table 1: NFW and Moore et al. density profiles without and with adiabatic compression 
(NFW C and Moore c respectively) with the corresponding parameters, and values of J(AO). 
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Figure 2: Dominant ncutralino annihilation diagrams. Relevant parts of the amplitudes are shown explicitly 
Terms between parenthesis correspond to /„ and Z final states in second and fourth diagrams. V and Z are 
chargino and neutralino mixing matrices. 



3 Neutralino dark matter 

For the computation of the gamma-ray flux discussed in the previous section, in addi- 
tion to the halo profile it is also crucial the value of the annihilation cross section cij. 
Of course, for determining this value the theoretical framework must be established. 

As mentioned in the Introduction, we will work in the context of the MSSM. Let us 
recall that in this framework there are four neutralinos, Xi (* = 1 ; 2, 3, 4), since they are 
the physical superpositions of the fermionic partners of the neutral electroweak gauge 
bosons, called bino (B°) and wino (W® ), and of the fermionic partners of the neutral 
Higgs bosons, called Higgsinos (H®, H®). Thus one can express the lightest neutralino 
as 

X? = ^n^ + Z 12 W$ + Z 13 H° d + Z lA H Q u . (3.7) 

It is commonly defined that x? is mostly gaugino-like if P = \Zn\ 2 + \Zu\ 2 > 0.9, 
Higgsino-like if P < 0.1, and mixed otherwise. 

In figure 2 we show the relevant Feynman diagrams contributing to neutralino 
annihilation. As was remarked e.g. in Ref. [26, 50, 51, 21], the annihilation cross section 
can be significantly enhanced depending on the SUSY model under consideration. We 
will concentrate here on the SUGRA scenario, where the soft terms are determined 
at the unification scale, Mqut ~ 2 x 10 16 GeV, after SUSY breaking, and radiative 
electroweak symmetry breaking is imposed. 

3.1 Supergravity models 

Let us discuss first the mSUGRA scenario, where the soft terms of the MSSM are 
assumed to be universal at Mqut- Recall that in mSUGRA one has only four free 
parameters: the soft scalar mass m, the soft gaugino mass M, the soft trilinear coupling 
A, and the ratio of the Higgs vacuum expectation values, tan/3 = (H®)/(Hj). In 
addition, the sign of the Higgsino mass parameter, (i, remains also undetermined by 
the minimization of the Higgs potential. 

Since in this scenario the lightest neutralino x? is mainly bino, only Z\\ is large 
and then the contribution of diagrams in Fig. 2 will be generically small, the first being 
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suppressed by / mass. As a consequence, the predicted annihilation cross section <jj 
is small, and therefore the flux is generically below the present accessible experimental 
regions (unless tan (5 > 50 combined with a Moore et al. profile are considered) 

However, as discussed in detail in Ref. [26] in the context of indirect detection, o; L 
can be increased in different ways when the structure of mSUGRA for the soft terms is 
abandoned. In particular, it is possible to enhance the annihilation channels involving 
exchange of the CP-odd Higgs, A, by reducing the Higgs mass. In addition, it is also 
possible to increase the Higgsino components of the lightest neutralino .£13,14. Thus an- 
nihilation channels through Higgs exchange become more important than in mSUGRA. 
This is also the case for Z~ , X\i and Xi-exchange channels. As a consequence, the 
gamma-ray flux will be increased. 

In particular, the most important effects are produced by the non-universality of 
Higgs and gaugino masses. These can be parameterised, at Mqut, as follows 

m 2 Hd =m 2 (l + 5 1 ) , m 2 Hu = m 2 (l + 5 2 ) , (3.8) 

and 

Mi = M , M 2 = M(l + S 2 ) , M 3 = M(l + 5' 3 ) , (3.9) 
We will concentrate in our analysis on the following representative cases: 
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Case a) corresponds to mSUGRA with universal soft terms, cases b), c) and d) cor- 
respond to non-universal Higgs masses, and finally cases e) and f) to non-universal 
gaugino masses. The cases b), c), d), and e) were discussed in Ref. [26], and are known 
to produce gamma-ray fluxes larger than in mSUGRA, whereas case f) will be of inter- 
est when discussing heavy WIMP signals predictions in the perspective of atmospheric 
Cherenkov telescopes like e.g. CANGAROO. 

We will discuss all these cases taking into account in the computation of the gamma- 
ray flux the effect of the adiabatic compression. Except for the calculation of the J 
factor, for the evaluation of the fluxes we used the last DarkSusy released version 
[52]. To solve the renormalization group equations for the soft SUSY-breaking terms 
between Mqut and the electroweak scale, we used the Fortran package SUSPECT [53]. 

Depending on the non-universal case, the main parameters that enter in the flux 
computation, namely the neutralino mass m x and the annihilation rate N^(av) = 

YIiIe^ dEy^ L (aiv) can vary significantly within the parameter space. Running a 
scan on the gaugino mass parameter M between and 2 TeV, we expect a neutralino 
mass lying between 200 GeV and 1 TeV, the lower bound being mainly restricted by 
accelerator constraints on the chargino or lighter higgs mass. We can however reach 
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3 to 4 TeV neutralinos respecting all the experimental and cosmological constraints if 
we extend the range of the scanning on M to 10 TeV. On the other hand, the WMAP 
relic density constraint restricts the annihilation cross section ^(av) to lie from ~ 
1 x 10 _28 cm 3 s _1 (when the coannihilation process with the lighter stau dominates, 
cases a) or c) with low tan/3) to ~ 2 x 10~ 26 cm 3 s _1 (when the A-pole region dominates, 
case c), d) and e) with tan/3 = 35) 3 . 

Before carrying out the analysis in detail, it is worth remarking that in these studies 
it is crucial to reproduce the correct phenomenology. Thus, in the next section, the most 
recent experimental and astrophysical constraints which can affect the computation will 
be discussed. 

3.2 Experimental and astrophysical constraints 

We have taken into account in the computation the most recent experimental and 
astrophysical bounds. These may produce important constraints on the parameter 
space of SUGRA models, restricting therefore the regions where the gamma-ray flux 
has to be analyzed. 

In particular, concerning the astrophysical bounds, the effect of 0.1 < Q,DMh 2 < 0.3, 
on the relic neutralino density computation has been considered. Due to its relevance, 
the WMAP narrow range 0.094 < QjjMh 2 < 0.129 has also been analyzed in detail. 
Actually, we have also considered the possibility that not all the dark matter is made 
of neutralinos, allowing 0.03 < Q^oh 2 < 0.094. In this case we have to rescale the 

Al 

density of neutralinos in the galaxy p(r) in Eq. (2.6) by a factor f2^o/i 2 /0.094. 

We illustrate this effect in Fig. 3 for the case of mSUGRA, i.e. case a) of (3.10). 
There, the cosmologically allowed regions are shown in the parameter space (m, M) for 
tan (3 = 50, A = and fi > 0. As discussed above, in most of the parameter space the 
lightest neutralino is mainly bino, and as a consequence the annihilation cross section 
is small producing a too large relic abundance. Nevertheless, let us recall that there 
are basically four corridors where the above bounds can be satisfied. There is the 
narrow coannihilation branch of the parameter space, i.e. the region where the stau is 
the next to the LSP producing efficient coannihilations (for 4 500 < m^o < 800 GeV). 
This reduces considerably the relic abundance and place it inside the bounds. Point 
A in the figure corresponds to this situation. Nevertheless, the neutralino annihilation 
in the galactic center is dominated by the A— Higgs exchange, but being far from the 
A— pole region (2m^o ~ mjCj the gamma-ray flux is not very large compared with the 
ones produced in the other regions to be discussed below (See Fig. 6a below, where 
this narrow corridor can be clearly distinguished). 

Now, for a given value of m, and from large to small M, itia decreases. When itia be- 
comes close to 2m ; o the annihilation cross section increases and therefore the relic den- 

Xi 

sity decreases entering in an allowed corridor (see e.g. point B) with 0.094 < ^^oh 2 < 

3 In some fine tuned regions of parameter space with especially large values of m (Focus Point region) it 
is possible to reach {av) ~ 10~ 25 cm 3 s _1 in the cases b), d) and e) 

4 Recall that we can deduce these values in the plot from the value of M, since in mSUGRA m^o ~ Mi ~ 
0.4 M. 



11 




500 1000 1500 2000 

M (GeV) 

Figure 3: Experimental and astrophysical bounds in the parameter space of the mSUGRA scenario (m, 
M) for tan/? = 50, A = and // > 0. The region to the left of the light grey (green) dashed line is excluded 
by the lower bound on the Higgs mass. The region to the left of the dotted line is excluded by the lower 
bound on the chargino mass m~± > 103 GeV. The region to the right of the black dashed line corresponds 
to a^ USY < 7.1 x 10~ 10 , and would be excluded by e + e~ data. The region to the left of the solid line is 
excluded by b — > 57. The region at the bottom (Stau LSP) is excluded because the lightest stau is the 
LSP. The region between grey (red) contours fulfils 0.1 < fi^o/i 2 < 0.3, whereas that between black contours 
indicates the WMAP range 0.094 < Q^oh 2 < 0.129. 
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0.3. Large gamma-ray fluxes dominated by bb final states through A— exchange are 
produced. For slightly smaller values of M one finally arrives to the A— pole region 
where rriA ~ 2m^o and the relic density is too small. Smaller values of uia close this 
region, producing a decrease in the annihilation cross section, which allows entry into 
the second corridor with the relic density inside the observational bounds (see e.g. 
point C). For iua too small this second corridor is closed, and the relic density is too 
large, ^^oh 2 > 0.3. 

Finally, there is the narrow Higgsino (focus-point) corridor (see e.g. point D) close 
to the no electroweak symmetry breaking (nEWSB) region where fi 2 becomes negative. 
In particular, the value of the relic density Q^o is affected due to the increase of the 
Higgsino components of x? with respect to the dominant bino component of most of the 
parameter space. Thus the relic abundance is placed inside the astrophysical bounds. 

The neutralino being mainly Higgsino couples strongly to the Z boson (see Fig. 2) 
giving rise to a large gamma-ray flux. For these points the neutralino is lighter than 
in the previous three cases. A full analysis of neutralino annihilation in the universal 
case parameter space can be found in [54]. 

Let us finally remark that qualitatively similar allowed corridors are found for the 
non- universal cases of (3.10). A detailed discussion can be found in Ref. [26]. 

Concerning the experimental constraints, the lower bounds on the masses of SUSY 
particles and on the lightest Higgs have been implemented, as well as the experimental 
bounds on the branching ratio of the b — > process and on a^ USY . We also illustrate 
these constraints in Fig. 3. Note that we are using \i > 0. We will not consider in the 
calculation the opposite sign of fi because this would produce a negative contribution 
for the — 2, and, as will be discussed below, we are mainly interested in positive 
contributions. Recall that the sign of the contribution is basically given by //M2, and 
that M, and therefore M2, can always be made positive after performing an U(1)r 
rotation. 

For a^ USY , we have taken into account the recent experimental result for the muon 
anomalous magnetic moment [55], as well as the most recent theoretical evaluations of 
the Standard Model contributions [56]. It is found that when e + e~ data are used the 
experimental excess in (<7„ — 2) would constrain a possible SUSY contribution to be 
G SUSY = (27.l±10)xl0~ 10 . At 2a level this implies 7.1 x 10~ 10 < a^ USY < 47.1 xl0~ 10 . 
It is worth noticing here that when tau data are used a smaller discrepancy with the 
experimental measurement is found. In order not to exclude the latter possibility we 
will discuss the relevant value a^ USY = 7.1 x 10~ 10 throughout the paper. For example, 
in Fig. 3 the region of the parameter space to the right of this (black dashed) line will 
be forbidden (allowed) if one consider electron (tau) data. 

On the other hand, the measurements of B — ► X s ^f decays at CLEO [57] and BELLE 
[58], lead to bounds on the branching ratio b — > 57. In particular we impose on our 
computation 2.33 x 10 -4 < BR(b — > S7) < 4.15 x 10 -4 , where the evaluation is carried 
out using the routine provided by the program micrOMEGAs [59]. This program is 
also used for our evaluation of a^ USY and relic neutralino density. 

Finally, we have considered in the computation the experimental limit on the B s — ► 
branching ratio [60]. Although the upper bound on this process B(B S — > /U + /i~) 
< 2.9 x 10~ 7 does not constrain the parameter space of mSUGRA, it has been stressed 
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Figure 4: Gamma-ray flux for a threshold of 1 GeV versus B(B S — > /i+^~) in case aj mSUGRA of Eq. (3.10) 
in the text, and non universal Higgs cases b) 5\ = 0,(52 = lj and <5i = —1,^2 = 0. Points depicted with 
light grey (red) points have 0.129 < Q^oh 2 < 0.3, with black points have 0.094 < Q^ali 2 < 0.129, and finally 
with dark grey (blue) points have 0.03 < fl^oh 2 < 0.094 with the appropriate rescaling of the density of 
neutralinos in the galaxy as discussed in the text. 
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recently that in the non-universal cases giving rise to large direct detection cross sec- 
tions the constraints can be very important [39, 40]. The reason being that there is a 
strong correlation between both observables. The branching ratio as well as the cross 
section increase for low values of the Higgs masses (and large tan 0) . 

However, in the case of indirect detection of neutralinos through gamma-ray fluxes 
the correlation is more diluted. In particular, when the neutralino has an important 
Higgsino component, the Z— exchange channel dominates and as a consequence there 
is no a direct link between the value of the gamma-ray flux and B(B S — ► ^ + ^~). We 
can see this in Fig. 4, where the gamma-ray flux versus the branching ratio is shown 
for cases a), b), c) of (3.10). As mentioned above, for mSUGRA the parameter space 
is not constrained. Although this is not the case for the non-universal examples, note 
that still large values of the gamma-ray fluxes can be obtained fulfilling the upper 
bound for B(B S — > 

4 Comparison with the experiments 

In this section we will compare the theoretical predictions for the indirect detection 
of neutralinos through gamma-ray fluxes, with the sensitivity of present and future 
experiments. As discussed in the Introduction, several experiments have already data 
that apparently cannot be explained with the usual gamma-ray background. Given this 
situation the comparison between the theoretical computation and the experimental 
data is crucial. To carry this out we will apply the formulas of Sect. 2.1 to the halo 
models with baryonic compression discussed in Sect. 2.2, in the context of a general 
SUGRA theory where non-universal soft terms can be present. We will discuss first the 
case of space-based detectors such as EGRET and GLAST, and finally atmospheric 
telescopes such as CANGAROO and HESS. 

4.1 EGRET 

The EGRET telescope on board of the Compton Gamma-Ray Observatory has carried 
out the first all-sky survey in high-energy gamma-rays (« 30 MeV - 30 GeV) over a 
period of 5 years, from April 1991 until September 1996. As a result of this survey, it has 
detected a signal [31] above about 1 GeV, with a value for the flux of about 10~ 8 cm~ 2 
s _1 , that apparently cannot be explained with the usual gamma-ray background 5 . The 
source, possibly diffuse rather than pointlike, is located within the 1.5° (Afl ~ 10~ 3 
sr) of the galactic center. Due to the lack of precision data in the high energy bins, 
it seems impossible however to distinguish any annihilation channel leading to this 
photon excess. 

In recent papers [15, 21] studies of the perspective of interpreting the spectral 
features of EGRET data as produced by neutralino annihilation in the galactic center 6 , 
were carried out. In particular, in Ref. [15] an interesting analysis was carried out in 

5 Although alternative explanations for this result have been proposed modifying conveniently the standard 
theory of galactic gamma ray [61], these turn out to be also controversial [62]. 
6 see also [29] for a general analysis. 
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the context of a toy dark matter model assuming an annihilation cross section av ~ 
3 x 10 -26 cm 3 s -1 scaling with the inverse of the relic abundance. Also neutralino dark 
matter in the mSUGRA context was considered, including experimental constraints 
coming from accelerator physics. In Ref. [21] the analysis was carried out in the 
framework of an effective MSSM with non-universal gaugino masses, obtaining that 
the EGRET excess cannot be reproduced with a NFW profile. 

We will show in this section that it is possible to reproduce the spectral features 
of EGRET data in the framework of general SUGRA, when a NFW profile including 
adiabatic compression, NFW C , is considered. 

The results can be seen in Fig. 5, where several cases of Eq. (3.10) have been 
studied for a scan on m and M from to 2 TeV and three values of tan/3 (5, 35 and 50). 
Only points of the parameter space producing interesting gamma-ray fluxes concerning 
the EGRET experimental data are shown. Concerning accelerator constraints, these 
points fulfil the lower bounds on the masses of SUSY particles and Higgs boson, as 
well as the experimental bounds on the branching ratios of the b — > and B s — > 
processes. They also fulfil WMAP bounds discussed in the previous section. In 
particular, the thin solid (blue) line corresponds to a NFW C density profile. Concerning 
the background of the diffuse gamma-ray flux of the inner galaxy at the energy range 
of interest for our analysis (from 100 MeV to 10 GeV), the main production of gamma- 
rays is the interaction of cosmic rays (mainly protons and helium nuclei) with the 
interstellar medium (atomic and molecular hydrogen, and helium). Neutral pions 7r° 
produced in this process radiate gamma-rays with a spectrum peaked at ~ 70 MeV 
(tott/2), dropping at high energies with an energy power law which follows the initial 
cosmic ray spectrum of the form E~ a with a ~ 2.7 [15, 62]. The normalization 
factor has been fixed to have the best fit in agreement to the standard scenario, fa = 
2 x 10- 6 cm- 2 s- 1 GeV -1 sr- 1 for 2 GeV < E hin < 4 GeV [15]. 

As one can see, with small values of tan/3 we are not able to find points reproducing 
EGRET data. However, interesting points can be found for larger values of tan/3. 
In particular, we show in the figure the cases tan/3=35,50 and one can see that with 
non-universal cases EGRET data can be reproduced. This is obtained for neutralino 
masses between 150 and 600 GeV. 

Let us remark that it is possible to differentiate each point of the parameter space 
(m, M) by its gamma-ray spectrum. Indeed, just a look at Fig. 5b for tan/3 = 35 shows 
two kind of spectra. The two highest lines correspond to points with (m, M) = (533 
GeV, 400 GeV) and (733 GeV, 533 GeV), and they are located inside the Higgsino 
region dominated by // final states through Z— exchange (see e.g. point D of Fig. 
3). On the other hand, the lower line corresponds to (m, M) = (266 GeV, 666 GeV), 
which lies inside the stau coannihilation region (see e.g. point A of Fig. 3), where the 
annihilation is weak but dominated by A-exchange, and the spectrum is thus a bb one. 
In the case of Fig. 5c with tan/3 = 35, the higher fluxes corresponds to the closing of 
the A— pole (see e.g. point C of Fig. 3) whereas the lower flux spectrum is obtained 
through the opening of the A— pole (see e.g. point B of Fig. 3). The same remarks 
can be made for the other cases illustrated in Fig. 5. 

We have also carried out the same analysis for a Moore et al. profile including 
adiabatic compression, Moore c , as shown in the figure with a dashed line for the differ- 
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Figure 5: Gamma-ray spectra <^ 1 {E 1 ) from the galactic center as functions of the photon energy for the 
SUGRA cases discussed in Eq. (3.10) for tan/3 = 5,35,50, A — and /i > 0, compared with data from the 
EGRET experiment. NFW and Moore et al. profiles with adiabatic compression are used with Ail ~ 10~ 3 
sr. All points shown after a scan on m and M (0 - 2000 GeV) fulfil the accelerator constraints discussed in 
the text, and WMAP bounds. 
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ent points of the parameter space. Clearly, due to the very singular behaviour of this 
profile many points are constrained by EGRET data. 

4.2 GLAST 

Crucial information on the origin and nature of the EGRET excess at the galactic 
center will be provided by another space-based detector, GLAST. This is scheduled 
for launch in 2007, and will perform an all-sky survey covering a larger energy range 
(« 10 - 300 GeV) than EGRET, with a wider effective area and a better resolution 
in energy. The experiment will also allow to search for an angular distribution of the 
dark matter halo. Indeed, with an angular resolution of 0.1° (AO ~ 10~ 5 sr), GLAST 
will be able to point and analyse the inner center of the Milky Way (~ 7 pc). Even if 
alternative explanations solve the data problem pointed out by EGRET, GLAST will 
be able to detect a much smaller flux of gamma rays from dark matter ~ 10~ n cm -2 

We summarize the results of our analysis in Fig. 6, using the same parameter 
space as in Fig. 5, and a NFW C profile. There, the values of the total gamma-ray 
flux $ 7 for a threshold of 1 GeV allowed by all experimental constraints as a function 
of the neutralino mass m^o are depicted. Concerning the predicted sensitivity for 
the experiment, we follow Ref. [65] and show on the plot the 5<r sensitivity curves 
assuming an angular resolution of 0.1° (AO = 10~ 5 ), and a total effective pointing 
time of 30 days for the satellite. We observe that, basically, the whole parameter 
space of general SUGRA will be tested by GLAST. Note that even for the mSUGRA 
case points corresponding to a value of tan/3 as small as 5 will be reached by GLAST. 
This is a remarkable result if we realise that interesting direct detection experiments, 
such as e.g. EDELWEIS-II or CDMS-II, will only be able to cover small regions of the 
parameter space. More poweful detectors such as GENIUS will cover larger regions, but 
still not too large compared with the parameter space of general SUGRA. In contrast, 
we have observed that indirect detection experiments will be able to test the parameter 
space of SUGRA 7 . 

It is also worth noticing that many points are already constrained by EGRET data. 
See e.g. case c) with tan (5 = 35. Let us remark that in Fig. 6 EGRET line corresponds 
to AO ~ 10~ 3 sr, and therefore for its study all points in the figure should be rescaled 
by a factor J(10~ 3 sr)10~ 3 / J(10 -5 sr)10~ 5 ~ 1.5. Of course, a Moore c profile will give 
rise to larger fluxes. 

4.3 CANGAROO 

Recently, the CANGAROO-II atmospheric Cherenkov telescope has made a signif- 
icant detection of gamma rays from the galactic center region [32]. In particular, 
the collaboration has published the spectrum obtained in six energy bins, from 200 
GeV to 3 TeV. Observations taken during 2001 and 2002 have detected a statisti- 
cally significant excess at energies greater than ~ 250 GeV, with an integrated flux 

7 See also Ref. [28] for a comparison between direct and indirect dark matter search without using 
pressed halo models. 
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Figure 6: Scatter plot of the gamma-ray flux $ 7 for a threshold of 1 GeV as a function of the neutralino 
mass m x for the SUGRA cases discussed in Eq. (3.10), and for the same parameter space used in Fig. 5. A 
NFW profile with adiabatic compression is used with AS1 ~ 1CP 5 sr. All points in the figure fulfil the b — > 57 
bounds. For tan /3=5 all points have a^ USY < 7.1 x 10~ 10 . For the other values of tan/3, a^ USY > 7.1 x 10~ 10 
is fulfilled for the points located on the arrow side of the dashed lines, and B(B S — > /i + /i~) < 2.9 x 10 -10 
for those on the arrow side of the dot-dashed lines. Points depicted with light grey (magenta) triangles have 
0.129 < %o/i 2 < 0.3, with black stars have 0.094 < Q^oh 2 < 0.129, and finally with dark grey (blue) boxes 
have 0.03 < ^^oh 2 < 0.094 with the appropriate rescaling of the density of neutralinos in the galaxy as 
discussed in the text. Solid lines denote the 5cr sensitivity curves for satellites. From top to bottom, the first 
solid line corresponds to the signal reported by EGRET, and the upper area bounded by the second solid 
line will be analyzed by GLAST experiment. 




Figure 7: Left frame: Gamma-ray spectra from the galactic center as functions of the photon energy for 
the case f) discussed in Eq. (3.10) with non-universal gaugino masses Mi = 0.5M, for tan/3 = 35, A = 
and /i > 0, compared with data from the CANGAROO-II experiment. NFW and Moore et al. profiles, as 
well as the average profile defined in Sect. 2.2 with adiabatic compression are used with Aft ~ 10~ 5 sr. All 
points shown after a scan on m and M (0-10 TeV) fulfil the accelerator constraints discussed in the text, 
and WMAP bounds. Right frame: Experimental and astrophysical bounds in the parameter space (m, M). 
The conventions and constraints are the same as in Fig. 3. Blue lines are contours of neutralino masses. 



of ~ 2 x 10~ 10 photons cm~ 2 s . These measurements indicate a very soft spectrum 
oc £)- 4 - 6±0 - 5 . As discussed in Ref. [63], this signal could be fitted with a spectrum of a 
TeV-dark matter candidate. It is interesting to see whether it is possible to obtain such 
a candidate in SUGRA scenarios imposing the accelerator and WMAP constraints, and 
withing the framework of adiabatically compressed halos. 

As we can see in Fig. 3, most of the regions of the parameter space of mSUGRA 
which satisfy WMAP constraints give mainly a sub- TeV spectrum in scalar and gaugino 
sectors. It is thus very difficult to find a dark matter candidate fitting the observations 
made by the CANGAROO-II experiment. Nevertheless, it is possible to extend the 
allowed region if we release the gaugino universality. One way to achieve this is to 
decrease M2 while keeping M\ = M3 at Mqut- This allows to reduce significantly 
the relic abundance of the wino-like neutralino, through its coannihilation with the 
lightest chargino X\ ■ Such an scenario, well motivated e.g. in anomaly mediation 
models, predicts higher values for the gaugino masses fulfilling WMAP constraints 
than those of universal cases, typically around the TeV scale, which is precisely the 
region favoured by CANGAROO-II results. 

As an illustration, let us consider the case f) in Eq. (3.10) giving rise to M2 = 0.5M. 
We show this case in Fig. 7. There all points fulfil the lower bounds on the masses of 
SUSY particles and Higgs boson, as well as the experimental bounds on the branching 
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ratios of the b — > S7 and -B s — > fi + processes. They also fulfil the WMAP bounds. 
As we can see in the left and right frames of the figure, there is an interesting region of 
the parameter space fulfilling the WMAP constraints and giving rise to a sufficiently 
large flux to fit the data from CANGAROO-II. Let us remark that the highest energy 
bin shown (of about 2.5 TeV) is less than 1 a in excess of a null result, so it should 
be considered only as an upper limit. Excluding this point from the analysis, we see 
that this non-universal model can fit CANGAROO data and WMAP bounds in halos 
models with profiles between NFW C and Moore c . In particular, this is the case of the 
scenario discussed in Sect. 2.2, NFW£., with a value of J given by 5 x Jnfw c > as can 
be seen in the left frame of the figure. 

It is worth noticing that the parameter space where WMAP constraints are fulfilled 
corresponds to a narrow range of the neutralino mass (~ 2 TeV), as can be seen in the 
right frame of Fig. 7. This region is mainly independent of tan/3 and the value of 5' 2 as 
far as the neutralino is wino like (5' 2 < —0.45). In this region of the parameter space, 
the relic density calculation is governed by the Xi ~ Xi coannihilation into W + W~ 
{m o ~ m + ~ M 2 ), whereas the flux is largely dominated by the annihilation process 

X? X? — * W~ through t— chanel xf exchange (cn;|wiy ~ 10 _26 cm 2 s _1 ). 

Another possibility is to produce a neutralino-LSP with a higher higgsino com- 
ponent. Indeed, as discussed in Ref. [26], decreasing the value of M3 while keeping 
M 2 = Mi at Mqut increases the Higgsino fraction of the neutralino through the 
renormalization group equations (RGEs), and, as a consequence, its coupling to the Z 
boson. In this case, the relic density constraint can be satisfied more easily because the 
contribution from the annihilation into gauge bosons channel is more important and 
the focus point region is much wider. Another consequence of the decreasing of M3 at 
Mqut is the decreasing of the pseudoscalar mass through m# 2 . In the left frame of 
Fig. 8 we show the case e) of Eq. (3.10) with M3 = 0.5M. For tan/3 larger than about 
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Figure 9: The same as in the left frame of Fig. 7 but for the non-universal case e) with M3 = 0.5M, 
compared with data for EGRET and CANGAROO-II experiments and the expected GLAST sensitivity. 
Here only the average profile defined in Sect. 2.2 with adiabatic compression, NFW^., is used. 



10, the neutralino annihilation into fermions (mainly bb) starts to dominate. Eventu- 
ally, the focus point and Higgs annihilation regions merge, for example for tan/3 = 35 
and M ~ 2 TeV, as was already pointed out by the authors of Ref. [64]. This region is 
clearly visible in the right frame of Fig. 8. Whereas in the previous analised case f) the 
WMAP constraint was fulfilled through one main chanel (xi Xi coannihilation), the 
non-universality of M3 offers more possibilities through the exchange of Higgses. That 
explains why the region allowed by WMAP constraint is broader in the right frame of 
Fig. 8 than in Fig. 7. 

Of course in these models the SUSY spectrum is very heavy. For example, for a 
point with m = 5 TeV, M = 5 TeV reproducing CANGAROO-II data, the spectrum in 
case e) (f) consists of heavy squarks, rriq ~ 8 TeV (10 TeV), neutralinos and charginos 
m o ~ m + ~ 1 TeV (2 TeV), and also heavy higgses itla ~ 5 TeV (6 TeV). If we do not 
take into account the potential fine-tuning problem associated with such models, and 
consider them reliable, this would be a difficult situation for any hope of discovering 
SUSY at colliders, since the only detectable particle would be the lightest Higgs. We 
have plotted for information the 2a standard deviation a^ USY = 7.1 x 10~ 10 with 
dashed line in the right frames of Figs. 7 and 8. Only an unambiguous evidence for a 
non-zero contribution to 8a ^ would restrict this class of models with heavy sparticles. 

Let us finally mention that we have also try to reproduce the flux measured by 
CANGAROO-II using cases b), c) and d) with non universal scalars. They give rise to 
a neutralino too light to be able to reproduce those data. 

4.4 Combining EGRET and CANGAROO 

After analysing EGRET and CANGAROO-II data, it seems natural to try to fit both 
experiments with only one non-universal scenario. We notice here that a similar ex- 
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ercise has been carried out by the authors of Ref . [65] . The baryonic cooling effect on 
the fluxes gives us the order of magnitude needed to fit with both data with a 1 TeV 
neutralino in the non-universal case e) with M3 = 0.5M. It is worth noticing that 
the CANGAROO-II collaboration in [32] pointed out already that the EGRET and 
CANGAROO-II data can be relatively smoothly connected with a cutoff energy of 1-3 
TeV. In our case, for that we need to use the NFW' C scenario discussed in Sect. 2.2. 
Typical points of the parameter space fullfilling all experimental constraints and fitting 
both set of data lie between (m = 800 GeV, M = 800 GeV) and (m = 3 TeV, M = 3 
TeV). 

It is also interesting to see the complementarity of GLAST with EGRET and CAN- 
GAROO. GLAST will perform an all-sky survey detection of fluxes with energy from 
1 GeV to 300 GeV, exactly filling the actual lack of experimental data in this energy 
range (see Fig. 9), and checking the CANGAROO results. Indeed, we have calculated 
that the integrated gamma ray flux for such a signal will be around 5 x 10~ n cm~ 2 s _1 . 
We have shown this sensitivity curve in Fig. 9 for AS7 = 10 -5 , which is the typical de- 
tector acceptance, following the prescriptions given in [65]. We clearly see that GLAST 
will help to cover the entire spectrum. 

We have also tried to reproduce the observed data using a wino-like neutralino aris- 
ing from the non-universal case f) with M2 = 0.5M. But any model of this type fulfill- 
ing WMAP bounds gives rise to a too heavy neutralino-LSP (~ 2 TeV) whose spectral 
features cannot explain at the same time both the CANGAROO-II and EGRET excess. 



4.5 HESS 

The HESS Cherenkov telescope experiment has recently published new data on gamma 
rays, detecting a signal from the galactic center [33]. The measured flux and spectrum 
differ substantially from previous results, in particular those reported by the CAN- 
GAROO collaboration, exhibiting a much harder power-law energy spectrum with 
spectral index of about —2.2 and extended up to 9 TeV. The authors of [33] already 
pointed out that if we assume that the observed gamma rays represent a continuum 
annihilation spectrum, we expect m x > 12 TeV. Actually such a heavy neutralino-LSP 
is not natural in the framework of a consistent supergravity model when we impose 
the renormalisation group equations and radiative electroweak symmetry breaking. 

Although we performed some scans in all non universality directions using numerical 
dichotomy methods, no point in the parameter space in any non-universal case studied 
was able to give a several 10 TeV neutralino satisfying WMAP constraint but this can 
be sensitive to the RGE and relic density calculation codes. 

Nevertheless, without RGE and taking soft parameters at the electroweak scale, 
the constraints are easier to evade. On top of that, in a very effective approach using 
completely free parameters and couplings in cross sections, neutralinos with m x > 
10 TeV and £l x h 2 ~ WMAP may certainly be fine tuned. We did not adopt such 
approaches since the MSSM is motivated by high energy and theoretical considerations. 
In this section, we analyze in a quite model-independant way the conditions required on 
the particle physics field to fit with the HESS data thanks to dark matter annihilation. 
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Figure 10: PYTHIA simulations for gamma spectrum coming from the decays of dark matter 
particle annihilation products ( a) bb, b) ti, c) W + W~ , d) t + t~ ) for three dark matter particle 
masses: 1,15 and 30 TeV. We show the data, our fit results and also compare with previous studies 
[13, 65]. 24 



We performed a simulation with the PYTHIA 8 package [66] to evaluate the gamma- 
ray fluxes from hadronisation and radiative processes. We consider a dark matter 
particle annihilating at rest in the bb, ti, W + W~ and t + t - final states. We first 
require a 1 TeV mass (i.e. E cm =2 TeV) to check our results with previous studies 
[13, 65] and then generate two relevant masses for HESS data: 15 TeV and 30 TeV 
(i.e. respectively E cm =30 and 60 TeV). The results are shown on figure 10. Even if 
more precise fits can be obtained with more parameters, the annihilation spectra can 
be quite simply approximated by : 

where x = E 1 /m x and we obtained (a, b) = (1.2, 10.) for bb 9 , (a, b) = (1.27, 14) for tt 
and (a, b) = (0.82,8.25) for W + W~ . For r + r~ we used the function proposed in [65] 

= x a ^ (b T x + c T x 2 + d T x 3 )e~ eTX (4.12) 

ax 

to well reproduce the harder r + r~ spectrum and found (a T ,b T ,c T ,d T ,e T ) = (—1.15,10,—' 
2,5). 

The gamma ray flux was expressed in Eq. (2.5), with (aiv) = BRi(av), BRi being 
the branching ratio of annihilation channel i. 
With this parameterization, the gamma flux can be rewritten as 



*7,_-».-.„..,^ _ 0.94 X10-". / („) \ (lOOOjVV 

m x (GeV) ^ dx \ 10- 29 cm 3 S - 1 ) \ m x J 



dE 1 



(4.13) 

where m x , (av), BRi are now considered as free parameters and the astrophysics 
contribution J is fixed by our halo density profiles. The results are shown in Fig. 11. 

We point the fact that we do not only compare the spectrum shape of the signal 
with possible dark matter annihilation explanation. Indeed, it should be noticed that 
our compressed halo profiles give rise to absolute gamma fluxes within the HESS data 
order of magnitude with (av) values in agreement with the WMAP requirement. 

According to our analysis, interpreting the HESS data in terms of dark matter 
particle annihilation is possible for a mass in the range 10 TeV < m x < 30 TeV and 
annihilation should be dominated by bb, ti or W + W~ final states with a possible 
t + t~ contribution 10 . The values for (av) required for NFW compressed and Moore 
compressed are respectively around 10 -25 cm 3 s~ 1 and 10~ 27 cm 3 s _1 . 

We also tried to see if it was possible to explain at the same time both HESS 
and EGRET excess by invoking a very heavy neutralino. We can easily predict the 
differential flux around 1 GeV for a 15 TeV neutralino. Using Eq. (4.11) for a it distri- 
bution (the one producing the less energetic photons) normalised with the HESS data 



8 We used january 2005 release of PYTHIA (version 6.317). 

9 Despite the Hill function, the bb case can suggest to relaxe the -1.5 exponent to get better fit. 
10 In the MSSM, ~ \^ ~ 0J5 for P scudo scalar exchange. 
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Figure 11: Dark matter annihilation versus HESS data for : NFW compressed in a) and b), Moore 
et al. compressed in c) and d). 
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in the 2.6 TeV energy bin (2. x l(r 16 cm- 2 s^GeV- 1 ) we obtain : d/V*/dE 7 (lGeV) ~ 
7.1 x 10 cm" - 2 s _1 GeV _1 , which is far below the EGRET excess data in the 1 GeV 
energy bin (1.52 x 10~ 7 cm~ 2 s~ 1 GeV _1 ). At the same time, we have numerically 
checked that such a heavy neutralino cannot be invoked to explain the EGRET excess 
in any of its energy bin. During this work, we saw that in a recent note added in [65], 
the authors reach the same conclusion. 

5 Conclusions 

We have analysed the effect of the compression of the dark matter due to the infall 
of baryons to the galactic center on the gamma-ray flux. In addition, we have also 
consider the effect of non-universal soft terms. This analysis shows that neutralino dark 
matter annihilation can give rise to signals largely reachable by future experiments like 
GLAST. In particular, we observe that, basically, the whole parameter space of general 
SUGRA will be tested by GLAST. Even for the mSUGRA case, points corresponding to 
a value of tan/3 as small as 5 will be tested. This is a remarkable result if we realise that 
direct detection experiments will only be able to cover a small region of the parameter 
space. In contrast, we have observed that indirect detection experiments will be able 
to test the parameter space of SUGRA. Let us point out that in this analysis we have 
taken into account all the recent experimental constraints, such as the lower bound on 
the Higgs mass, the b — ► 57 branching ratio, the muon g — 2 and the recently improved 
upper bound on B(B S — > / Lt + ^~), as well as the astrophysical (WMAP) bounds on the 
dark matter density. 

Actually, in this SUGRA framework we have also been able to fit present excess 
from EGRET and CANGAROO using different non-universal scenarios, and even fit 
the data from both experiments with only one scenario. We have also studied the 
recent HESS data implying a neutralino heavier than 12 TeV. Because of such a heavy 
neutralino, it is not natural to find solutions in the SUGRA framework. Nevertheless we 
have carried out a quite model-independent analysis, and found the conditions required 
on the particle physics side to fit the HESS data thanks to dark matter annihilation. 

In any case, we must keep in mind that the current data obtained by the different 
gamma-rays observations from the galactic center region by these three experiments 
do not allow us to conclude about a dark matter annihilation origin rather than other, 
less exotic astrophysics sources. Fortunately, this situation may change with the im- 
provement of angular resolution and energy sensitivity of future detectors like GLAST, 
which will be a complementary source of gamma-ray data. 
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